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Abstract. The main proposal in this paper is the merging of two tech- 
niques that have been recently developed. On the one hand, we consider 
a new approach for computing some specializable Grobner basis, the so 
called Minimal Canonical Comprehensive Grobner Systems (MCCGS) 
that is -roughly speaking- a computational procedure yielding "good" 
bases for ideals of polynomials over a field, depending on several param- 
eters, that specialize "well", for instance, regarding the number of solu- 
tions for the given ideal, for different values of the parameters. The sec- 
ond ingredient is related to automatic theorem discovery in elementary 
geometry. Automatic discovery aims to obtain complementary (equality 
and inequality type) hypotheses for a (generally false) geometric state- 
ment to become true. The paper shows how to use MCCGS for automatic 
discovering of theorems and gives relevant examples. 
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1 Introduction 

1.1 Overview of Goals 

The main idea in this paper is that of merging two recent techniques. On the one 
hand, we will consider a method (named MCCGS, standing for minimal canon- 
ical comprehensive Grobner systems) [MaMoOGj, that is -roughly speaking- a 
computational approach yielding "good" bases for ideals of polynomials over a 
field depending on several parameters, where "good" means that the obtained 
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bases should specialize (and specialize "well" , for instance, regarding the number 
of solutions for the given ideal) for different values of the parameters. 

Briefly, in order to understand what kind of problem MCCGS addresses, let 
us consider the ideal (ax,x + y)K[a][x, y], where a is taken as a parameter and 
if is a field. Then it is clear that there will be, for different values of a = ao G K, 
essentially two different types of bases for the specialized ideal (ciqx, x+y)K[x, y]. 
In fact, for ao = we will get (x + y) as a Grobner-basis (in short, a G-basis) 
for the specialized ideal; and for any other rational value of a such that a = 
a 7^ 0, we will get a G-basis with two elements, (x, y). Thus, the given G-basis 
{ax, x + y)K[a, x, y] does not specialize well to a G-basis of every specialized 
ideal. On the other hand, let us consider (ax — b)K[a, b][x], where a, b are taken 
as free parameters and x is the only variable. Then, no matter which rational 
values ao, b a are assigned to a, 6, it happens that {agx — b } remains a Grobner 
basis for the ideal (a x — b )K[x]. Still, there is a need for a case-distinction if 
we focus on the cardinal of the solutions for the specialized ideal. Namely, for 
a ^ there is a unique solution x = —b /a ; for a = and b ^ there is 
no solution at all; and for ao = bo = a solution can be any value of x (no 
restriction, one degree of freedom). 

The goal of MCCGS is to describe, in a compact and canonical form, the dis- 
cussion, depending on the different values of the parameters specializing a given 
parametric system, of the different basis for the resulting specialized systems 
and on their solutions. 

The second ingredient of our contribution is about automatic theorem discov- 
ery in elementary geometry. Automatic discovery aims to obtain complementary 
hypotheses for a (generally false) geometric statement to become true. For in- 
stance, we can consider an arbitrary triangle and the feet on each sides of the 
three altitudes. These three feet give us another triangle, and now we want to 
conclude that such triangle is equilateral. This is generally false, but, under what 
extra hypotheses (of equality type) on the given triangle will it become generally 
true? 

Finding, in an automatic way, the necessary and sufficient conditions for this 
statement to become a theorem, is the task of automatic discovery. A proto- 
col for automatic discovery is presented in [RV99] and a detailed discussion of 
the method appears in |DR| . The protocol proceeds requiring some computa- 
tions (contraction, saturation, etc.) about certain ideals built up from the given 
statement, but does not state any preference about how to perform such compu- 
tations (although the computed examples in both papers rely on straightforward 
Grobner bases computations for ideal elimination). 

Our goal in this paper is to show how we can improve the automatic discovery 
of geometry theorems, by performing a MCCGS procedure on an ideal built 
up from the given hypotheses and theses, considering as parameters the free 
coordinates of some elements of the geometric setting, 




1.2 Related work 

This idea has a close precedent in the worljf] of |CLLW| . directly inspired by 
[K95J and, to a lesser extent, by [Wcis92J. In [CLLWj . a parametric radical 
membership test is presented for a mathematical construct the authors intro- 
duce, called "partitioned parametric Grobner basis" (PPGB). Suppose we are 
given a statement H := {hi = 0, . . . , h r = 0} T := {g = 0}, expressed in 
terms of polynomial equations -usually over some computable field- and their 
solutions over some extension held K -that we can assume, in order to simplify 
the exposition, to be algebraically closed. Roughly speaking, the method behind 
[CLLWJ starts by computing the "partitioned basis" (with respect to a given 
subset of variables, here denoted by u) of an ideal I C K [u, x, y], (for instance, 
/ = (/ii(u,x) . . . h r (u,x),g(u, x) y — 1)), ie. a finite collection of couples (Ci,Fi), 
where the C^'s are constructible sets described as {ci = 0, . . . , c m = 0,qi ^ 
0, . . . , q s 0} on the parameter space, and the Fi's are some collections of poly- 
nomials in K[u, x, y]. Moreover, it is required (among other conditions) that the 
Cj's conform a partition of the parameter space and, also, that for every element 
u in each d, the (reduced) G-basis of (/ii(u ,x) . . . h r , (u , x), g(u , x) y — 1) 
is precisely Fi(u ,x). It is well known (e.g. |K86j or |Ch88j ) that, in this con- 
text, a statement {hi = . . . h r = 0} {g = 0} is to be considered true if 
1 e (hi . . . h r , gy — 1); thus, the extra hypotheses that |CLLW| proposes to add 
for the statement to become a theorem are precisely those expressed by any of 
the Ci's such that the corresponding Fi — {1}, since this is the only case Fi can 
specialize to {1}. 

We must remark that, simply testing for 1 G (h\ . . . h r , g y — 1), as in the 
method above, can yield to theorems that hold just because the hypotheses are 
not consistent (i.e. such that already 1 S (hi . . . h r ) ). This cannot happen with 



But notice the authors of [CLLWj already mention the paper of Montes |Mo02| as a 
predecessor on this particular kind of discussion of Grobner basis with parameters. 



our approach to automatic discovery: if a new statement is discovered, then the 
obtained hypotheses will be necessarily consistent. 

Although our approach stems from the same basic ideas, our contribution 
differs from [CLLWj in some respects: first, we focus on automatic discovery, 
and not in automatic proving. Moreover, we are able to specifically describe the 
capability and limitations of the method (while in |CLLW| it is only mentioned 
that, in the reducible case, their "method . . . cannot determine if the conclusion 
of the geometric statement is true on some components of the hypotheses" ) . Sec- 
ond, even for proving, the use of MCCGS provides not only the specialization 
property (which is the key for the application of partitioned parametric bases 
in [CLLWj ) but also an automatic case distinction, that allows a richer under- 
standing of the underlying geometry for the considered situation. In fact, it 
seems that the partitioned parametric G-Basis (PPGB) algorithm from [CLLW 
is close to the algorithm DISPGB considered in [Mo02j, both sharing that their 
output requires collecting by hand multiple cases (and then having to manually 
express in some simplified way the union of the corresponding conditions on 
the parameters). Actually the motivation for MCCGS was, precisely, improving 
DISPGB. 

Our approach has also an evident connection (since |Weis92| is the common 
origin of all posterior developments on parametric Grobner basis) to the work 
of several members of Prof. Weispfenning's group, regarding generic quantifier 
elimination (Q. E.) and its application to automatic theorem proving (as, for 
example, in |DG| . [DSWj . [SSj . [St]). In particular we remark the strong relation 
of our work with that of [DG| , that approaches theorem proving via a restricted 
(generically valid) Q.E. method, relying on generic Grobner systems computa- 
tions. The set of restrictions provided by this method, besides speeding up 
the Q.E. computations, can be interpreted in the context of theorem proving, 
roughly speaking, as a collection of new (sufficient) non-degeneracy conditions 
for an statement to hold true. 

Again, the difference between our contribution here and theirs is, first, that 
we address problems requiring, in general, parameter restrictions that go beyond 
"a conjunction & of negated equations in the parameters" ( jDGj . first paragraph 
in Section [3]) . That is, we deal with formulas that are almost always false (see 
below for a more detailed explanation of the difference between automatic deriva- 
tion and automatic discovery) and require non-negated (ie. equality) parameter 
restrictions; they can not be directly approached via generic Q.E. since our for- 
mulas are, quite often, generically false. Moreover, our approach is limited to 
this specific kind of generically false problems and we do not intend to provide 
a general method for Q.E. A second difference is that, for our very particular 
kind of problems, MCCGS formulates parameter restrictions in a compact and 
canonical way, a goal that is not specifically intended concerning the description 
of in |DG] . For these reasons we can not include performance comparisons to 
these Q.E. methods and we do not consider relevant (although we provide some 
basic information) giving hardware details, computing times, etc. on the perfor- 
mance of our method running on the examples described in the last section of 



this paper. We are not proposing something better, but something different in a 
different context. 

Next Section includes a short introduction to the basics on automatic discov- 
ery, which could be of interest even for automatic proving practitioners. Section 
[3] provides some bibliographic references for the problem of the G-basis special- 
ization and summarizes the main features of the MCCGS algorithm, including 
an example of its output. Section 0] describes the application of MCCGS to auto- 
matic discovery, while Section 5 works in detail a collection of curious examples, 
including the solution of a pastime from Le Monde and the simpler solution (via 
this new method) of one example also solved by a more traditional method. 

2 A digest on automatic discovery 

Although less popular than automatic proving, automatic discovery of elemen- 
tary geometry theorems is not new. It can be traced back to the work of Chou 
(see |Ch84j . |Ch87| and |ChG90] l. regarding the "automatic derivation of for- 
mulas" , a particular variant of automatic discovery where the goal consists in 
deriving results that always occur under some given hypotheses but that can be 
formulated in terms of some specific set of variables (such as expressing the area 
of a triangle in terms of the lengths of its sides). Finding the geometric locus of 
a point defined through some geometric constraints (say, finding the locus of a 
point when its projection on the three sides of a given triangle form a triangle of 
given constant area |Ch88j . Example 5.8) can be considered as another variant 
of this "automatic derivation" approach. 

Although "automatic derivation" (or locus finding) aims to discover some new 
geometric statements (without modifying the given hypotheses) , it is not exactly 
the same as "automatic discovery" (in the sense we have presented it in the 
previous section), that searches for complementary hypotheses for a (generally 
false) geometric statement to become true (such as stating that the three feet 
of the altitudes for a given triangle form an equilateral triangle and finding 
what kind of triangles verify it). Again, automatic discovery in this precise sense 
appears in the early work of Chou (whose thesis |Ch85j deals with "Proving and 
discovering theorems in elementary geometries using Wu's method") and Kapur 
[K89] (where it is explicitly stated that ". . . the objective here is to find the 
missing hypotheses so that a given conclusion follows from a given incomplete 
set of hypotheses. . . "). 

Further specific contributions to automatic discovery appear in |Wa98j , [R98 
(a book written in Spanish for secondary education teachers, with circa one 
hundred pages devoted to this topic and with many worked out examples), 
[RV99], [Koj or |CWj . Examples of automatic derivation, locus finding and dis- 
covery, achieved through a specific software named GDI (the initials of Geometria 
Dindmica Inteligente), of Botana-Valcarce, appear in |BR05| or |RB| (and the 
references thereof), such as the automatic derivation of the thesis for the cele- 
brated Maclane 83-Theorem, or the automatic answer to some items on a test 



posed by Richard [Ri , on proof strategies in mathematics courses, for students 
14-16 years old. 

The simple idea behind the different approaches i^jj, essentially, that of adding 
the conjectural theses to the collection of hypotheses, and then deriving, from 
this new ideal of theses plus hypotheses, some new constraints in terms of the 
free parameters ruling the geometric situation. For a toy example, consider that 
x — a — is the only hypothesis, that the set of points (x, a) in this hypothesis 
variety is determined by the value of the parameter a, and that x = is the 
(generally false) thesis. Then we add the thesis to the hypothesis, getting the new 
ideal (x — a,x), and we observe that the elimination of x in this ideal yields the 
constraint a — 0, which is indeed the extra hypothesis we have to add to the given 
one x — a = 0, in order to have a correct statement [x — a = A a = 0] => [x = 0] . 

With this simple idea as starting pomt0, an elaborated discovery procedure, 
with several non trivial examples, is presented in [RV99J. It has been recently 
revised in [BDR] and |DR] , showing that, in some precise sense, the idea of 
considering H + T for discovering is intrinsically unique (see Section [4] for a 
short introduction, leading to the use of MCCGS in this context). 

3 Overwiew on the MCCGS algorithm 

As mentioned in the introduction, specializing the basis of an ideal with param- 
eters does not yield, in general, a basis of the specialized ideal. 

This phenomenon -in the context of Grobner basis- has been known for over 
fifteen years now, yielding to a rich variety of attempts towards a solution (we 
refer the interested reader to the bibliographic references in |MaMo06] or in 
|Wib06j). Finding a specializable basis (ie. providing a single basis that collects 
all possible bases, together with the corresponding relations among the parame- 
ters) is -more or less- the task of the different comprehensive G-Basis proposals. 
Although the first global solution was that of Weispfenning, as early as 1992 
(see [Wcis92 ), the topic is quite active nowadays, as exemplified in the above 
quoted recent papers. The MCCGS procedure, that is, computing the minimal 

4 Already present in the well known book of Ch88], page 72: "...The method 
developed here can be modified for the purpose of finding new geometry theo- 
rems. . . Suppose that we are trying to prove a theorem. . . and the final remain- 
der. . . Ro is nonzero. If we add a new hypotheses Ro — 0, then we have a theo- 
rem. . . " . Here Chou proposes adding as new hypotheses the pseudoremainder of the 
thesis by the ideal of hypotheses, a mathematical object which should be zero if the 
theorem was generally true. 

5 Indeed, things are not so trivial. Consider, for instance, H => T, where H = (a + 
l)(a + 2)(6 + 1) C K[a, b, c] and T = (o + b + 1, c) C K[a, b, c]. Take as parameters 
U — {b, c}, a set of dim(H )-variables, independent over H. Then the elimination 
of the remaining variables over H + T yields H' = (c, b 3 — b). But H + H' — 
(a + 1, b, c) n (a + 2, b, c) H (a + 1, b - 1, c) n (o + 2, 6- 1, c) n (6 + 1, c) does not imply 
T, even if we add some non-degeneracy conditions expressed in terms of the free 
parameters U, since T vanishes over some components, such as (a + 2, b — 1, c) (and 
does not vanish over some other ones, such as (o + 1, b — 1, c)). 



canonical comprehensive Grobner system of a given parametric ideal, is one of 
the approaches we are interested in. Let us describe briefly the goals and output 
of the MCCGS algorithm. 

Given a parametric polynomial system of equations over some computable 
field, such as the rational numbers, our interest focuses on discussing the type of 
solutions over some algebraically closed extension, such as the complex numbers, 
depending on the values of the parameters. Let x = (x\, . . . ,x n ) be the set of 
variables, u = (ui, . . . , u m ) the set of parameters and / C K[u] [x] the parametric 
ideal we want to discuss, where, in order to simplify the exposition, a single field 
K, algebraically closed, is considered both for the coefficients and the solutions. 
We want to study how the solutions over K n of the equation system defined 
by / vary when we specialize the values of the parameters u to concrete values 
Uo € K. Denote by A = K[u], and by a Uo : A[x] — > K[x\ the homomorphism 
corresponding to the specialization (substitution of u by some u € K). 

A Grobner System GS(I, y x ) of the ideal / C A[x] wrt (with respect to) 
the termorder ^ x is a set of pairs (Si,Bi), where each couple consists of a 
constructible set (called segment) and of a collection of polynomials, such that 



The algorithm MCCGS (Minimal Canonical Comprehensive Grobner Sys- 
tem) |Mo06j . M;iMol)(i of the ideal / c A[x\ wrt the monomial order ^ x for the 
variables, builds up the unique Grobner System having the following properties: 

1. The segments Si form a partition S — {S\, . . . , S s } of the parameter space 
K m . 

2. The polynomials in Bi are normalized to have content 1 wrt x over K[u] 
(in order to work with polynomials instead of with rational functions). The 
Bi specialize to the reduced Grobner basis of <r Uo (J), keeping the same Ipp 
(leading power products set) for each u £ Si, i.e. the leading coefficients 
are different from zero on every point of Si . 0. Thus a concrete set of Ipp can 
be associated to a given Si- Often it exists a unique segment corresponding 
to each particular Ipp, although in some cases several such segments can 
occur. In any case, when a segment with the reduced basis [1] exists, then it 
is unique. When two segments Si, Sj share the same Jpp, then there is not a 
common reduced basis B specializing to both Bi, Bjj. 




Moreover, there exists a unique segment S± (called the generic segment), con- 
taining a Zariski-open set, whose associated basis B\ is called the generic 
basis and coincides with the Grobner basis of I considered in J^(u)[x] con- 
veniently normalized without denominators and content 1 wrt x. 

6 The polynomials in the B^s are not faithful (they do not belong to I), as they are 
reduced wrt to the null conditions in Si. By abuse of language we call them reduced 
bases (i.e. not-faithful, in the terminology of Weispfenning) . 

7 M. Wibmer Wib06] has proved that for homogeneous ideals in the projective space 
there is at most a unique reduced basis and segment corresponding to a given Ipp. 



GS(J, >- x ) - {(S u Bi) : 1 < i < s, Si C K m , B t C A[x\, (J* Si = K m 
Vuo S Si, a Uo {Bi) is a Grobner basis of a Uo (I) wrt ^ x }- 




3. The partition S is canonical (unique for a given / and monomial order). 

4. The partition is minimal, in the sense it does not exists another partition 
having property 2 with less sets Si. 

5. The segments Si are described in a canonical form. 

As it is known, the Ipp of the reduced Grobner basis of an ideal determine 
the cardinal or dimension of the solution set over an algebraically closed field. 
This makes the MCCGS algorithm very useful for applications as it identifies 
canonically the different kind of solutions for every value of the parameters. This 
is particularly suitable for automatic theorem proving and automatic theorem 
predicting, as we will show in the following sections. 

Let us give an example of the output of MCCGS. 

Example 1. Consider the system described by the following parametric ideal 
(here the parameters are a, b, c, d): 

I = (x 2 + by 2 + 2cxy + 2dx, 2x + 2cy + 2c?, 2by + 2cx), 

arising in the context of finding all possible singular conies and their singular- 
ities. Calling to the Maple implementation of MCCGS yields a graphical and 
an algebraic output. The graphical output is shown in Figure It contains the 
basic information that is to be read as follows. At the root there is the given 
ideal (in red). The second level (also in red) contains the Ipp of the bases of 
the three different possible cases. These are [1], corresponding to the no solution 
(no singular points) case; [x, y], corresponding to the one solution (one singular 
point) case; and [x], corresponding to the case of one dimensional solution (ie. 
when the conic is a double line). Below each case there is a subtree (in blue) 
describing the corresponding Si, with the following conventions: 

— at the nodes there are ideals of ii"[u], prime in the field of definition (gen- 
erated over the prime field by the coefficients of a reduced G-Basis) of the 
given ideal I c A[x] 

— a descending edge means the set theoretic "difference" of the set defined by 
the node above minus the set defined at the node below, 

— nodes at the same level, hanging from a common node, are to be interpreted 
as yielding the set theoretic "union" of the corresponding sets; they form the 
irredundant prime decomposition of a radical ideal of -K"[u]. 

— every branch contains a strictly ascending chain of prime ideals. 

So, in the example above, the three cases, their Ipp and the corresponding 
Si's are to be read as shown in the following table: 



Ipp 


Basis Bi 


Description of Si 


[1] 


[1] 


K :i \ ((V(6) \ (V(c, 6) \ V(d, c, &))) U V(d)) 


[y,x\ 


[2cy + d,x] 


(V(fo)\V(c,6))u(V(d)\V(d,6-c 2 )) 


[x] 


[x + cy] 


Y(d,b-c') 



We remark that the Bi 's do not appear in the Figure [2 since -in order to 
simplify the display- the complete bases are only given by the algebraic output 
of MCCGS and are not shown by the graphic output. 




4 Using MCCGS for automatic theorem discovering 

Once we have briefly described the context for MCCGS and for automatic dis- 
covery, we are prepared to describe the basic idea in this paper. We can say 
that our goal is to show how performing a MCCGS procedure can improve the 
automatic discovery of geometry theorems. 

Example [1] can be seen as a very simple example of theorem discovering. We 
could formulate the statement a conic has one singular point and try to find the 
conditions for the statement to be true. Without loss of generality we express 
the equation of the conic and its partial derivatives as 

/ = 2 + by 2 + 2cxy + 2dx, 2x + 2cy + 2d, 2by + 2cx), 

and search for the values of the parameters where this system has a single so- 
lution. As shown above, we have found that the statement is true if and only if 
{b = 0, c 7^ 0} or if {d — 0, b — c 2 ^ 0}, since in the first segment of the table 
there is no solution (£?i = (1)), while the third segment yields a 1-dimensional 
set of solutions. 

In general, let H T be a statement expressed in terms of polynomial 
equations, where the ideals iJ, T C K[x\, . . . , x n ] will be the corresponding 
hypotheses ideal and theses ideal (both, possibly, with several generators). In 
this context [PR] sets, as discovery goal, finding a couple of subsets of variables 
U 2 U', with X D U 2 U', and a couple of ideals R' C K[U],R" C K[U'], 
so that the following properties hold for the associated algebraic varieties (over 
K n , with K algebraically closed): 

1. Y(H + R' e )\Y(R" e ) C V(T) (where the e stands for the extension of the 
ideal from its defining ring, say, K[U] or K[U'], to K[X]); 

2. V(H + T) C V(i?' e ); 

3. Y(H + R' e )\Y(R" e ) ^ 0. 



The rationale behind such a definition is that such a couple (i?', R") is sup- 
posed to provide 

— some necessary (as expressed by item 2) above) 

— and sufficient (as expressed by item 1) above) 

non trivial (as expressed by item 3) ) complementary conditions of equality kind 
(given by R' ) and of non-degeneracy type (given by the negation of R") for the 
given theses to hold under the given hypotheses. 

Then it is shown in [PR] that, for a given couple of subsets of variables 
UDU', with X DU DU', there is a couple of ideals R' C K[U],R" C K[U'] 
verifying properties 1), 2) and 3) above if and only if the couple of ideals H' = 
(H + T) n K[U] ancO H" = ((H + H' e ) : T°°) n K[U'} also verify these three 
conditions. Moreover, Theorem 2 in |DR| shows these conditions hold if and 
only if 1 g {H') c : H' ,oc (equivalently, iff H" % ^{H') c ), where c stands for the 
contraction ideal, so there is an algorithmic way of solving the posed discovery 
problem for a given statement and choice of variables. 

Now we remark the following: 

Proposition 1. If there is a couple R' , R" verifying the above conditions, then 

Y(H + R /e )\Y(R" e ) = Y(H + H' e )\Y(R" e ) 

Proof. First notice that Y( H + H' e ) C Y(H + R' e ), sinc e, by pr operty 2), Y(H + 
T) C Y{R' e ), thus R' e C y/H + T and so R 1 = R' ec C y/H + T c = \fW and this 
implies that Y{H' e ) C Y{R' e ). 

Moreover, we have also that Y(H + R' e )\Y(R" e ) C V(T) n Y(H), by prop- 
erty 1), and V(T) n Y(H) C Y(H' e ), where the last inclusion follows from the 
definition of H' e . We conclude that Y(H + R' e )\Y(R" e ) C V(i? + H' e )\Y(R" e ). 

This means that the search for candidates i?' for complementary hypotheses 
of equality type, can be reduced to computing Y(H'). This is, precisely, the 
(Zariski closure of the) projection, over the parameter space of the U- variables, 
of Y(H+T), and this can be computed through MCCGS, providing as well some 
other useful information (as in Corollary [T]). 

Proposition 2. The projection ofY(H + T) over the U -variables can be com- 
puted by performing a MCCGS for I = H +T and X D U, discarding, if it exists, 
the unique segment Si 's with Bi equal to 1 and keeping the remaining Si 's. 

Proof. Since the segments of a MCCGS partition the parameter space U, it is 
enough to show that a point (uq) is not in the projection if and only if it belongs 
to the Si with associated Bi = 1 . Now we recall that a reduced Grobner basis is 
1 if and only if the corresponding ideal is (1). Then, the ideal H + T specialized 
at wo will be (1) if and only if its reduced G-basis is 1. Since we work over an 
algebraically closed field, this is the only case the system H+T, specialized at uq, 

8 Let I, J be ideals of -K'pf]. Recall that I : J = {x,xj C /}. Then, the saturation of 
I by J is denned as I : J°° = U„(J : J"), cf. [KROO] . 



has no solution, ie. uq is not in the projection of Y(H +T). But, by construction, 
a Bi specializes to 1 if and only if Bi — 1 (since the specialization must be a 
reduced G-basis and has the same lpp as Bi). 

Corollary 1. The union of these Si's with associated Bi ^ 1 (ie. the com- 
plement of the only possible segment with Bi = 1) partitions the projection of 
W(H + T); that is, it holds H AT =>• {US^}. Thus, the union of these Si 's provide 
complementary necessary conditions for the theses T to hold over H . 

We will see below (Remark [2]) that, when the given statement does not hold 
over any geometrically meaningful component of the hypotheses variety - ie. in 
the automatic discovery situation- the segment with Bi = 1 is the generic one, 
so its complement provides necessary conditions for the theses T to hold over 
H. 

Next we must study if some of these Si's provide sufficient conditions, ana- 
lyzing the behavior of each statement H A Si =>• T, for every segment Si with 
lpp =/= 1. Some -perhaps all, perhaps none- of them could be true. Remark 
that, anyway H A Si ^ 0, since the associated basis in not 1. Remark, also, 
that MCCGS allows to obtain supplementary conditions Si of the more general 
form (not every constructible set is the difference of two closed sets of the form 
Y(H + i?' e )\V(i?" e ), as in the previous approach). 

There are some special easy cases, as shown in the next result. 

Corollary 2. For every segment Si such that the corresponding lpp of the as- 
sociated basis is, precisely, the collection of variables {x\, . . . , x n }, we have that 
Y(H) n Si C V(T), ie. H A Si ^ T holds, and Si provides sufficient conditions 
for T to hold over H . 

Proof. In fact, the condition on the associated lpp means that for every uq in 
iSj, the system H(uq,x) = 0,T(uq,x) = has a unique solution, and it belongs 
to V(T). Thus Y(H) n S* 4 C V(T). 

Otherwise, we should analyze, for each i with Si involved in the projection 
of Y(H + T), the validity of H A Si =>• T. This is a straightforward "automatic 
proving" step, and not of "automatic discovery" , since adding again T to the 
collection of hypotheses H A Si will not change the situation, as the projection 
of Y(H) n V(T) n Si equals the projection of Y(H) n Si, both being 

Yet, MCCGS can provide a method for checking the truth of such statement 
H A Si => T. As it is well known, we can reformulate the hypotheses H A Si as a 
collection of equality hypotheses H, since Si is constructible and, then, the union 
of intersections of closed and open sets (in the Zariski topology). And open sets 
can be expressed through equalities by means of saturation techniques (such as 
x ^ <^ xy — 1=0, etc.). So let us state the following propositions (adapting 
to the MCCGS context some results from [RV99 , [DR|) in all generality. 

Proposition 3. Let H => T be a statement and let U be a collection of variables 
independent for H . Then T vanishes identically on all the components of H where 
U remain independent if and only if, performing a MCCGS for {H, Tz — 1} with 
respect to U , the generic basis is 1. 



Proof. Notice the stated condition on the segments of the MCCGS is equivalent 
to the fact that the contraction (H, Tz — 1)C\K[U] 7^ (0). In fact, this contraction 
is zero if and only if the projection of Y(H + (Tz — 1)) contains an open set. 
And this is equivalent to the fact that the generic segment has lpp 7^ 1. 

Now, if {H,Tz - 1) n K[U] ^ (0), take some ^ g e (H,Tz — 1) PI K[U]. 
Remark that, by construction, gT = over Y(H). If T 7^ at some point over 
some component of Y(H), then g = over such component; so it cannot be a 
component where the U are independent, since g € -K"[E/]. 

Conversely, if T vanishes identically over all the independent components, 
then we can compute an element g € K[U] vanishing over the remaining compo- 
nents (because U' is dependent over them). So gT vanishes all over Y(H), and 
thus (H, Tz ~ 1) n K[U'\ ^ (0). 

Remark 1. In fact, as in [CLLWj . it is easy to show that the segment with 
associated lpp equal to 1 provides complementary sufficient conditions for H =>■ T 
to hold. In fact, for every Uo in such segment, V(iJ(uo, x), T(uo, x)z — 1) = 0, 
so V(-ff(uo,x) C T(uo,x). But it can happen there is no such segment. 

Proposition 4. Let H =>■ T be a statement and let U be a collection of vari- 
ables independent for H and of dimension equal to dira(H). Then T vanishes 
identically on some components of H where U remains independent if and only 
if performing a MCCGS for {H, T} with respect to U , the reduced basis of the 
generic segment is different from 1. 

Proof. As above, the stated condition on the segments of the MCCGS is equiv- 
alent to the fact that the contraction (H,T) n K[U'] = (0). 

Now, if T does not vanish identically over any component of Y(H) indepen- 
dent over U' , the projection of Y(H,T) over U will be a proper closed subset 
(since the dimension of the projection is less or equal than the dimension of the 
components of Y(H) contained in Y(T) , the maximum dimension of all compo- 
nents of Y(H) equals the maximum dimension of the independent components, 
and the dimension of the [/-space equals the maximum dimension of the com- 
ponents of Y(H)). This contradicts the assumption (H,T) n K[U] = (0), which 
implies the closure of the projection is the whole [/-space. 

Conversely, if T vanishes identically over some independent component (say, 
C) and (H,T) n K[U] ^ (0), then we can choose an element 7^ g £ (H,T) n 
K[U]. This element vanishes over any component of Y(H) where T vanishes, in 
particular over C, contradicting its independence over U. 

Remark 2. The last proposition can be also read in a different way: T does not 
vanish identically on any independent component of H if and only if the reduced 
basis of the generic segment is 1. 

Corollary 3. Let H T be a statement and let U be a collection of vari- 
ables independent for H and of dimension equal to dim(iJ). Then T vanishes 
identically on some components of H where U remains independent and also T 
does not vanish identically on some other components of H where U remains 
independent if and only if 



— performing a MCCGS for {H, Tz — 1} with respect to U , the generic segment 
dos not have reduced basis 1, and 

— performing a MCCGS for {H, T} with respect to U , the reduced basis of the 
generic segment is also different from 1. 

In conclusion, using MCCGS one can determine, for a given statement, 
whether it is generally true (over all independent components, using Proposi- 
tion [31 generally false (over all independent components, using Remark [5]), or 
partially true and false (using Corollary [3]). Let us call this last situation the 
"undecidable" case. 

In fact, unfortunately, in this circumstance it is not possible, using only 
data on the U variables, to determine the components of H where T vanishes 
identically. Consider H = b(b + 1) C K[a,b], T = (b) and take U = {a}. Here 
the projection ofY(H,Tz — 1) over the U- variables is the whole a- line, so does 
not have any segment with Ipp equal to 1, and we know the thesis does not hold 
over all independent components. Moreover the projection oiV(H + T) over the 
{/-space is again the whole a-line, so there is no segment with Ipp 1, and we 
can conclude T holds over a component, but there is no way of separating the 
component b — 0, by manipulating H, T in terms of polynomials in the variable 
a. 

This discussion applies to the situation described above, when considering 
statements H A Si => T, where segment Si belongs to a MCCGS for {H, T} 
with respect to a collection U of variables and has Ipp ^ 1. Let HH be the 
reformulation of H A Si in terms of equalities and let (if possible) U' C U be 
a new collection of variables, such that they are independent for HH and of 
dimension equal to dira(HH). 

Then, as remarked above, HH T will be true on the segment SSi of a 
MCCGS with respect to HH,Tz — 1, with Ipp 1. If it is is an open segment, 
then the statement H A Si => T will be generally true (over all the components 
independent over U'). If it is not open segment, but there is at least one such 
segment, the statement will hold true under the new restrictions. 

But if there is no segment at all with Ipp = 1, then and only then we are in the 
undecidable case. In fact, over all points in the t/'-projection of Y(HH, Tz— 1) we 
will have points of Y(H) not in V(T) (because all the segments will have Ipp ^ 1 
in the MCCGS for HH, Tz-l) and also points of V(H) and V(T) (since we are 
also in the projection of Si over U', and Si corresponds to a segment of Ipp =/= 1 
for a MCCGS with respect to H, T). 

In this case, since the projection over U' of V(HH, T) will be same as the 
projection of Y(HH) (both being equal to the projection of Si), it is of no use 
to go further with a new discovery procedure, computing a MCCSG for HH, T 
over U'. We know beforehand that all its segments will have Ipp ^ 1, since over 
any point in the projection of Si there will be always points on Y(HH) n V(T), 
confirming, again, that we are in the undecidable situation. 



5 Examples 



Let us see how this works in a collection of examples, where we have just de- 
tailed the discovery step (ie. computing the MCCGS of {H, T} with respect to 
a collection of maximal independent variables for H, and then collecting the 
potentially true statements H A Si =>■ T, where segment Si has lpp 7^ 1) in the 
procedure outlined in the previous Section. That is, we have not included here 
the formal automatic verification in each case that the newly found hypotheses 
actually lead to a true statement (the "proving step"). 

Example 2. (See also [DRj ) . Next, we will develop the above introduced notions 
considering a statement from [Ch88 (Example 91 in his book), suitably adapted 
to the discovery framework. The example here is taken from [DR . 

Let us consider as given data a circle and two diametral opposed points on it 
(say, take a circle centered at (1, 0) with radius 1, and let C = (0, 0), D — (2, 0) 
the two ends of a diameter), plus an arbitrary point A — (1*1,1x2). See Figure[3J 
Then trace a tangent from A to the circle and let E = (xi, x%) be the tangency 
point. Let F — (x^, X4) be the intersection of DE and CA. Then we claim that 
AE = AF. Moreover, in order to be able to define the lines DE, CA, we require, 
as hypotheses, that D ^ E (ie. u\ 7^ 2) and that C ^ A (ie. u\ 7^ or u 2 7^ 0). 




Fig. 3. Problem of Example [2] 



Now, using CoCbA [CNR99] and its package TP (for Theorem Proving), we 
translate the given situation as follows 



Alias TP := $contrib/thmproving; 



Use R: :=Q[x[l. .4] ,u[l. . 2] ] ; 

A: = [u[l] ,u[2]] ; 
E: = [x[l] ,x[2]] ; 
D:=[2,0] ; 
F: = [x[3] ,x[4]] ; 
C:=[0,0] ; 

Ipl:=TP. Perpendicular ([E, A] , [E, [1,0]]) ; 
Ip2:=TP.LenSquare([E, [1,0]] )-l; 
Ip3:=TP.Collinear( [0,0] ,A,F) ; 
Ip4:=TP.Collinear(D,E,F) ; 

H : =Saturat ion (Ideal (Ipl , Ip2 , Ip3 , Ip4) , Ideal (u [1] -2) * 

Ideal(u[l] , u[2])) ; 

T : =Ideal (TP . LenSquare ( [A , E] ) -TP . LenSquare ( [A , F] ) ) ; 

where T is the thesis and H describes the hypothesis ideal. Notice that Ipl 
expresses that the segments [E, A], [E, (1,0)] are perpendicular; Ip2 states that 
the square of the length of [E, (1, 0)] is 1 ( so Ipl, Ip2 imply E is the tangency 
point from A ); and the next two hypotheses express that the corresponding 
three points are collinear. The hypothesis ideal H is here constructed by using 
the saturation command, since it is a standard way of stating that the hypothesis 
variety is the (Zariski) closure of the set defined by all the conditions Ip[i] = 
0,i = 1 . . .4 minus the union {u[l] — 2} U {u[l] = 0, u[2] = 0}, as declared in 
the formulation of this example (but we refer to |DR) for a discussion on the 
two possible ways of introducing inequalities as hypotheses). Finally, the thesis 
expresses that the two segments [AE] , [AF] have equal non oriented length. 

Now, let us use in this, clearly false, statement the approach of |RV99] or 
DR to discovery. First we check that the statement H — > T is not algebraically 
true in any conceivable way. For instance, it turns that 

Saturation(H, Saturation (H, T) ) ; 
Ideal (i) 



and this computation shows that all possible non-degeneracy conditions (those 
polynomials p(u, x) that could be added to the hypotheses as conditions of the 
kind p(u, x) ^ 0) lie in the hypotheses ideal, yielding, therefore to an empty set 
of conditions of the kind p ^ Ap = 0. This implies, in particular, that the same 
negative result would be obtained if we restrict the computations to some subset 
of variables, since the thesis does not vanish on any irreducible component of 
the hypotheses variety. 



Thus we must switch on to the discovery protocol, checking before hand that 
u[l],u[2] actually is a (maximal) set of independent variables -the parameters- 
for our construction: 



Dim(R/H) ; 
2 



Elim([x[l] ,x[2] ,x[3] ,x[4]] ,H); 
Ideal (0) 



Then we add the thesis to the hypotheses ideal and we eliminate all variables 
except u[2] 

H' :=Elim([x[l] ,x[2] ,x[3] ,x[4]] ,H+T) ; 
H'; 

Ideal (-l/2u[l] ~5 - l/2u [1] ~3u [2] ~2 + u[l]~4) 



Factor(-l/2u[l] ~5 - l/2u [1] ~3u [2] ~2 + u[l]~4); 
[[u[l]"2 + u[2]~2 - 2u[l], 1], [u[l], 3], [-1/2, 1]] 



yielding as complementary hypotheses the conditions u[l] 2 + u[2] 2 — 2u[l] = 
V u[l] = that can be interpreted by saying that either point A lies on the 
given circle or (when u[l] = 0) triangle A(A, C, D) is rectangle at C. In the 
next step of the discovery procedure we consider as new hypotheses ideal the set 
H + H', which is of dimension 1 and where both u[2] or u[l] can be taken as 
independent variables ruling the new construction. 

Dim(R/(H+H')) ; 
1 



Elim([x[l] ,x[2] ,x[3] ,x[4] ,u[l]] ,H+H>) ; 
Ideal (0) 



Elim([x[l] ,x[2] ,x[3] ,x[4] ,u [2] ] ,H+H' ) ; 
Ideal (0) 

Choosing, for example, u[2] as relevant variable, we check -applying the usual 
automatic proving scheme- that the new statement H A H ' — > T is correct under 
the non-degeneracy condition u[2] ^ 0: 

H> ' :=Elim([x[l] ,x[2] ,x[3] ,x[4] ,u[l]] , Saturation(H+H' ,T)) ; 
H"; 

Ideal (u [2] ~3) 



Thus we have arrived to the following statement: Given a circle of radius 1 and 
centered at (1, 0), and a point A not in the X-axis and lying either on the Y axis 
or in the circle, it holds that the segments AE, AF (where E is a tangency point 
from A to the circle and F is the intersection of the lines passing by (2,0), E 
and A, (0, 0)) are of equal length. 

Let us now review Example [5] using MCCGS. As above, the hypotheses are 
the union of H := Hi U S, where H\ expresses the equality type constraints: 

Hi = [(xi - l){ui - xi) + x 2 {u 2 - x 2 ), {xi - l) 2 + x% - I, 
U1X4 — U2X3, X3X2 — X4X1 — 2x2 + 2x4] 

to which we have to add the saturation ideal expressing the inequality con- 
straints: 

S = [U1X4, — U2X3, X1U1 — Iti — X\ + X2U2, X4X2 — 2X2U2 — X3U1 + 2ui, 

X4X1 — 2xiu 2 + u 2 x 3 , X3X2 — 2xiu 2 + u 2 x 3 — 2x 2 + 2x 4 , 
X1X3 + X3U1 + 2X2U2 — 2xi ~ 2ui,x\ — 2xi + x\, 
X3U1 + 2X2U2U1 — 2u\xi + U2X3 — 2u\ — 2X2U2 + 2112X4, 
x\u\ + X4U2X3 + 2x\ — 4x3^1 — AU2X4 + Aui, 
uix\ — X1X2U2 — x\ + X2U2 + xi — ui, 

U2X3 + U2X4X3 + 2x\ — AU2X\ — Au 2 x\ + AU2X3}. 

The thesis is 

T = (Ui - Xif + (u 2 - X 2 ) 2 - (Ui - X3) 2 - (U2 - X 4 ) 2 - 

Calling now mccgs(_ff! U S U T, lex(a;i, X2, x%, X4), lex(M!, m 2 )) one obtains the 
following segments: 



Segment 


lpp 


Description of Si 


1 


[1] 


K 2 \ (Y{uj +u 2 2 - 2ui) U Y{ui)) 


2 


[x 2 4,Xz,X2,Xl] 


V(ui + u 2 - 2m) \ (V(ui - 2, u 2 ) U V(ui, u 2 )) 


3 


[x'(,X 3 ,X2,Xi} 


Y{ui)\{Y{ui,ui + l)\JY{ui,u 2 )) 


4 


[Xi, X3, X2,Xl] 


Y(ui,u 2 + 1) 


5 


[xj,X3,xl,Xi] 


V(ui-2,ua) 


6 


[x 2 4,xl,x 2 ,xi] 


V(u 2 ,ui) 



Segment Si states that point A(ui,u 2 ) must lie either in the Y-axis or on 
the circle, as a necessary condition in the parameter space u = (1x1,112) for the 
existence of solutions, in the hypotheses plus thesis variety, lying over u. This 
essentially agrees with the result obtained in [DRj . 

A detailed analysis of the remaining segments show a variety of formulas for 
determining the (sometimes not unique) values of points E(xi, x 2 ) and F(x3, X4) 
-verifying the theorem- over the corresponding parameter values. 



For completeness we give the different bases associated, in the different seg- 
ments, to the above ideal of thesis plus hypotheses 

Bi = [1] 

B 2 = [u 2 + x\ — 2u 2 x 4 , — u\x 4 + U2X3, t*2 — 2u 2 Ui + x 2 u\ + {—2u\ + 2ui)x 4 , 

U 2 U\ + X1U2 — 2U1X4] 

B 3 = [-2u 2 x 4 + x% x 3 , (u 2 + l)x 2 - a; 4 , {u 2 + l)xi - 112X4] 

Bi = [x4,X 3 ,X 2 ,Xi] 

B 5 = [#4, -4 + 2x 3 , xl,-2 + xi] 

B& = [x\, x%, —x 3 x 4 + 2x 2 — 2x 4 , 2xi] 

Example 3. Next we consider the problem^) described in Figure [4] Take a circle 
C with center at 0(0,0) and radius 1 and let us denote points A = (—1,0) and 
B = (0, 1). Let D be an arbitrary point with coordinates D = (1 + a, b) and let 
C = (1 + a,0) be another point in the X-axis, lying under point D. Then trace 
the line BC. Assume this line intersects the circle C at point P(x,y). 



£>'(!+ a, &') 




Fig. 4. Example [3] 



Consider now the, in general false, statement "the points A, P, D are aligned". 
We want to discover the conditions on the parameters a, b for the statement to 



We thankfully acknowledge here that this problem was suggested by a colleague, 
Manel Udina 



be true. The set of hypotheses plus thesis equations are very simple: 

HT = [x 2 + y 2 - 1, -x + 1 - y + a - ay, -2y + b + xb-ay] 

Take x, y as variables and a, b as parameters and call mccgs(HT, lex(x, y), lex(o, b)). 
The graphical output of the algorithm can be seen in Figure O and the algebraic 
description appears in the following table. 



lpp 


Basis Bi 


Description of Si 


[1] 


[1] 


(K 2 \ (V(a - 6) \ V(a - b, (b + l) a + 1))) 
U (K 2 \V{2 + a)) 
U (K 2 \Y(a-b + 2)) 




[x' + y'-l, 

x + {a + l)(y-l), 

b(x + l)- (a + 2)y] 


(V(o-6)\V(o-6, (6 + l) a + l)) 
U (V(2 + o)\V(6,2 + o)) 
U (V(o-6 + 2)\V(ft,2 + a)) 


[y 2 ,x] 


[y(y- + 


¥(6, 2 + a) 




[1] 



Fig. 5. Canonical tree for Example [3] 



As we see, the generic case has basis [1] showing that the statement is false 
in general. The interesting case corresponds, as it is usually expected, to the 
case with lpp = [x, y], providing a unique solution for P. The description of the 
parameter set associated to this basis gives the union of three different locally 
closed sets, namely V(o - 6) \ V(a - b, (b + l) 2 + 1), V(2 + a) \ V(6, 2 + a) and 
V(a — 6+2)\V(6, 2 + a), expressing complementary hypotheses for the statement 
to hold. 

The first set is (perhaps) the expected one, corresponding to the case a = b 
(except for the degenerate complex point (6, b) with (b + l) 2 + 1 = 0, without 



interest from the real point of view). Thus we can say that the statement holds 
if point C is equidistant from point D and point E. 

The second set yields a = — 2 and corresponds to the situation where point D 
is on the tangent to the circle trough the point (—1,0) (except for the degenerate 
case b = 0). In this case P = A and, obviously, A, P, D are aligned (even in the 
degenerate stated in the third segment, corresponding to the lpp = 

[y 2 M)- 

Finally, the third set gives the condition b = a + 2 and it is also interesting, 
since it corresponds to the case where the intersecting point of the line BC with 
the circle is taken to be B instead of P, and then point D' should be in the 
vertical of C and at distance D'C equal to distance EC plus two. 

Example 4- [Isosceles orthic triangle] 

In |DRj the conditions for the orthic triangle of a given triangle (that is, 
the triangle built up by the feet of the altitudes of the given triangle over each 
side) to the equilateral have been discovered. Next example aims to discover 
conditions for a given triangle in order to have an isosceles orthic triangle. 

Consider the triangle of Figure[T]with vertices A{— 1, 0), B(l, 0) and C(a, b), 
corresponding to a generic triangle having one side of length 2. Denote by 
Pi(a, 0), P 2 (x 2 ,y 2 ), ^3 (#3, 2/3) the feet of the altitudes of the given triangle, 
ie. the vertices of the orthic triangle. The equations defining these vertices are: 

H=(a-l)y 2 -b(x 2 -l)=Q,' 
{a-l)(x 2 + l) + by 2 = 0, 
(o + l)|ft-6(aj 3 + l) = 0, 
(a+l)(x 3 -l) + by 3 = 0, ^ 



Now let us add the condition P1P3 = P\P 2 . 

T=(x 3 - a) 2 + y 2 - (x 2 - a) 2 - y 2 = 0. 
Take x 2 ,x 3 ,y 2 , y 3 as variables and a, b as free parameters and call 

mccgs(ff U T, lex(x2, x 3l y 2 , y 3 ), lex(a, b)). 

The output has now four segments. The generic case, with lpp = [1], meaning 
that the orthic triangle is, in general, not isosceles; one interesting case with 
lpp — [2/3, 2/2, X3, x 2 \\ and two more cases we can call degenerate, with Ipp's 
[y 2 ,x 2 ,x 2 ] and [2/2, £3, £2] > respectively. For the interesting case we show the 
graphic output in Figure [5] Its basis is 

B 2 = [(a 2 + b 2 + 2a+ l)y 3 - 2ab - 2b, (a 2 + b 2 - 2a + l)y 2 + 2ab - 2b, 

(a 2 + b 2 + 2a + l)x 3 - a 2 + b 2 - 2a - 1, (a 2 + b 2 - 2a + \)x 2 + a 2 - b 2 - 2a 

Next table shows the description of the lpp and the Si's for the the four cases: 



[y3, y2, x3, x2] 




[b A 2+1,a] [b A 2+1,a] [b, a-1] [b, a+1] [b, a-1] [b, a+1] 

Fig. 6. Canonical tree branch for lpp = [y 3 , j/ 2 , 2:3, a; 2 ] in Example!?] 



lpp 


Description of S 1 , 


[1] 


X 2 \((V(a)\V(6 2 + l,a)) 

U (V(a 2 -6 2 -l)\V(6 2 + l,a)) 

U V(a 2 + 6 2 - 1)) 


[?/3, 2/2, »3, X 2 ] 


V(a)\V(6 2 + l,a) 

U V(a 2 + 6 2 - 1) \ (V(6, a - 1) U V(6, a + 1)) 

U (V(a 2 - b 2 - 1) \ (V(& 2 + 1, a) U V(6, a - 1) U V(6, a + 1)) 


[j/2, 33,2:2] 


V(6,a + 1) 


[2/2,^3,^2] 


V(6,a- 1) 



The description of the parameter set (over the reals) for which the theorem is 
potentially true and no degenerate can be phrased as follows: 



1) a = 

2) a 2 + b 2 = 1 except the points (1, 0) and (-1, 0) 

3) a 2 — b 2 = 1 except the points (1,0) and (—1,0) 

This set is represented in Figure [7] and corresponds to 

1) The given triangle is itself isosceles (a — 0); 

2) The given triangle is rectangular at vertex C (with vertices A(—l, 0), B(l, 0) 
and the vertex C(a, b) inscribed in the circle a 2 + b 2 = 1, 

3) The given triangle has vertices A{— 1, 0), B(1,0) and vertex C(a,b) lies on 
the hyperbola a 2 — b 2 = 1. 

Solution 1) is, perhaps, not surprising. Solution 2) corresponds to rectangular 
triangles for which the orthic triangle reduces to a line, that can be considered 
a degenerate isosceles triangle. But solution 3) is a nice novelty: it exists a one 
parameter family of non-isosceles triangles having isosceles orthic triangles. 

The remaining two cases in the MCCGS output with lpp = [2/2,^3,^2] 
and lpp = [2/2,2:3,2:2] represent degenerate triangles without geometric inter- 
est (namely C = A and C = B). 

Thus, after performing an automatic proving procedure for the new hypothe- 
ses, we can formulate the following theorem: 



Fig. 7. Solutions of Example |4] 



Theorem 1. Given a triangle with vertices A(— 1, 0), B(1,0) and C(a,b), its 
orthic triangle will be isosceles if and only if vertex C lies either on the line 
a = (and then the given triangle is itself isosceles) or in the circle a 2 + b 2 — 1 
(and then it is rectangular) or in the hyperbola a 2 — b 2 = 1. 

Example 5. [Skaters] 

Our final example is taken from the pastimes section of the French journal Le 
Monde, published on the printed edition of Jan. 8, 2007. This example is there 
attributed to E. Busser and G. Cohen. We think it is nice from Le Monde to 
include the proof of a theorem as a pastime. Actually, the statement to be proved 
was presented as arising from a more down-to-earth situation: two ice-skaters are 
moving forming two intersecting circles, at same speed and with the same sense 
of rotation. They both depart from one of the points of intersection of the two 
circles. Then the journal asked to show that the two skaters were always aligned 
with the other point of intersection (where some young lady, both skaters were 
interested at, was placed...). 

Let us translate this problem into a theorem discovering question, as follows. 

We will consider two circles with centers at P(a, 1) and Q(—b, 1) and radius 
r\ = a 2 + 1 and r\ = b 2 + 1, as shown in Figured! intersecting at points O(0, 0) 
and M(0,2). Consider generic points -the skaters- A(x\,yi) and B(x2,yi) on 
the respective circles. Point A will be parametrized by the oriented angle v = 
OP A and, correspondingly, point B will describe the oriented angle w = OQB. 
Therefore we can say that angle zero corresponds to the departing location of 
both skaters, namely, point O. 

We claim that, for whatever position of points A, B, the points A, M, B are 
aligned, which is obviously false in general. But we want to determine if there is 
a relation between the two oriented angles making this statement to hold true. 



Denote c v ,s v ,c w ,s w the cosine and sine of the angles v and w. It is easy to 
establish the basic hypotheses, using scalar products: 

Hi = [(an - a) 2 + ( Vl - l) 2 - a 2 - 1, (x 2 + b) 2 + (y 2 - l) 2 - b 2 - 1, 

a(si - a) + (j/i - 1) + (a 2 + l)c„, -6(0:2 + 6) + (y 2 - 1) + (1 + b 2 )c w 

Now, as the angles are to be taken oriented (because we assume the skaters 




Fig. 8. Skaters problem 



tare moving on the corresponding circle in the same sense), we need to add the 
vectorial products involving also the sine to determine exactly the angles and 
not only their cosines. So we add the hypotheses: 

H 2 = [a(yi - 1) - (a* - a) + (a 2 + l)s v , -b(y 2 - 1) - (x 2 + b) + (b 2 + l)s w ] 

The thesis is, clearly: 

T = x\y 2 - 2x\ - x 2 yi + 2x 2 . 
The radii of the circles are 

r\ = a 2 + 1 and r\ = b 2 + 1 
and for n 7^ and r 2 ^ we have 

c Vn = cos v a = cos OPM = — , s Vn = sin vo = sin OPM = — — ^— , 

a z + 1 " a z + 1 

- 6 2 - 1 2b 

c wo = cosw = cosOQM = -, s wa = s'mw = smOQM = — -. 

o^ + l b z + 1 



We want to take a, b and the angles v and w -in terms of the sines and cosines- 
as parameters. So we must introduce the constraints on the sine and cosine pa- 
rameters. Moreover, we notice there are also some obvious degenerate situations, 
namely n = 0, r 2 = and a + b = 0, corresponding to null radii or coincident 
circles, and we want to avoid them. 

Currently, MCCGS allows us to introduce all these constraints in order to 
discuss the parametric system. The call is now 

mccgs(i/i U H 2 U T, lex(ari, yi, x 2 , 2/2), lex(a, b, s v , c v , s w , c w ), 

null = [c 2 + s 2 v - 1, c 2 + s 2 - 1], notnull = {a 2 + 1, b 2 + 1, a + b}). 

including the constraints on the parameters and eluding degenerate situations 
as options for MCCGS. 

The result is that MCCGS outputs only 2 cases. The first one has basis [1], 
showing that, in general, there is no solution to our query. The second one has 
lpp = [y 2 , x 2l yi, xi] determining in a unique form the points A and B for the 
given values of the parameters. The associated basis is 

[2/2 +c w — bs w -l,x 2 - bc w - s w +b,yi+ c v + as v — l,xi + ac v - s v - a] 

with parameter conditions that can be expressed as the union of three irreducible 
varieties: 

V 2 = V(c 2 u + s 2 w - 1, c 2 + si - 1, s w + bc w - 6, bs w - c w - 1) 
V 3 = V(c 2 u + - 1, cl + si - 1, -s v + ac v - a, as v + c v + 1) 

The interpretation is easy: V\ corresponds to arbitrary a, b, w, plus the essen- 
tial condition v = w, which is the interesting case, stating that our conjecture 
requires (and it is easy to show that this condition is sufficient) that both skaters 
keep moving with the same angular speed. 

V 2 corresponds to s w = s Wo , c w = c Wo and a, b, v free, thus B — M and A 
can take any position. 

V3 is analogous to V 2 , and corresponds to placing A — M and B anywhere. 

So we can summarize the above discussion in the following 

Theorem 2. Given two non coincident circles of non-null radii and centers P 
and Q, intersecting at two points O and M, let us consider points A, B on 
each of the circles. Then the three points A, M, B are aligned if and only if the 
oriented angles OP A and OQB are equal or A or B or both coincide with M . 

6 Performances 

Although the principal advantage of MCCGS in relation to other CGS algo- 
rithms is the simplicity and properties of the output: the minimal number of 
segments and the characterization of the type of the solution depending on the 



values of the parameters, the computer implementation^! of the corresponding 
package, named dpgb release 7.0, in Maple 8 is relatively short time consuming. 
Moreover, we think that no other actual PCAD software will be able to obtain 
the accurate result obtained, for example, in example 13. We give here a table 
with the CPU time and number of segments for the examples of the paper. 



Example 


CPU time (sec.) 


Number of segments 


m 


1.9 


3 


m 


12.8 


6 


m 


0.98 


3 


® 


4.4 


4 


m 


129.4 


2 



The computations were done with a Pentium(R) 4 CPU at 3.40 Ghz and 
1.00 GB RAM. 

7 Conclusion 

We have briefly introduced the principles of automatic discovery and also the 
ideas -in the context of comprehensive Grobner basis- for discussing polynomial 
systems with parameters, via the new MCCGS algorithm. Then we have shown 
how natural is to merge both concepts, since the parameter discussion can be 
interpreted as yielding, in particular, the projection of the system solution set 
over the parameter space; and since the conditions for discovery can be obtained 
by the elimination of the dependent variables over the ideal of hypotheses and 
thesis. Moreover, we have also remarked how the approach through MCCGS 
provides new candidate complementary conditions of more general type and, in 
some particular instances (segments of the parameter space yielding to unique 
solution), quite common in our examples, an easy test for the sufficiency of these 
conditions. Finally, the use of MCCGS for automatic proving has been presented, 
as part of a formal discussion on the limitations of the discovery method. 

We have exemplified this approach through a collection of non-trivial ex- 
amples (performed by running the current Maple implementation of MCCGS, 
see ]MaMo06 , over a laptop, without special time - a few seconds- or memory 
requirements), showing that in all cases, the MCCGS output is very suitable 
to providing geometric insight, allowing the actual discovery of interesting and 
new? theorems (and pastimes!). 
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